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Abstract 

We study a simple map as a minimal model of excitable cells. The map has two fast 
variables which mimic the behavior of class I neurons, undergoing a sub-critical Hopf 
bifurcation. Adding a third slow variable allows the system to present bursts and 
other interesting biological behaviors. Bifurcation lines which locate the excitability 
region are obtained for different planes in parameter space. 
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1 Introduction 



Maps are dynamical systems with continuous state variables but discrete-time 
dynamics. They constitute low coniputational cost elements for large scale 
simulations of complex systems The proper use of coupled map lattices 

(CML) in neural modeling depends on the choice of a good basic element. 
Two- variable maps have been proposed by Aihara et. al. Chialvo and 
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Kinouchi and Tragtenberg j5|. Methodological considerations about the use of 
maps in computational neuroscience can be found in these references. 

A three-variable map with one slow and two fast variables has been recently 
introduced by Kuva et al. |6]. It is able to reproduce complex biological be- 
haviors, such as bursting, cardiac-like spikes, chaotic orbits and slow regular 
spikes ("singlet bursting"). Here we determine stability diagrams for this map, 
which may be thought of as a discrete-time analogous of the Hindmarsh-Rose 
neuron [7|. We focus exclusively on the lines at which the fixed point loses 
stability, which define the limit of the excitability region. We extend previ- 
ous results by obtaining diagrams beyond the adiabatic limit, where the slow 
variable has a much slower time scale than the other two j6|. 



2 The model 

The following three-dimensional non-linear map has been proposed as a min- 
imal excitable bursting cell model 0: 



Like in the Hindmarsh-Rose model |7|, the variable x represents the instan- 
taneous membrane potential, y is a recovery variable and z can be consid- 
ered a slow adaptive current if 6 and A are set to small values (<^ 1). Under 
these conditions, we can identify z{t) with e.g. the Ca^^-dependent current 
Iahp [1, 1^, since its time scale and general effect (long lasting hyperpolariza- 
ton and adaptation of firing rates) are very similar. I{t) represents an external 
current which can model the effects of a synapse, a stimulus-induced ionic flow 
in sensory cells, or an experimentally controlled clamped current. 

If we set z{t) = H = const, we are left with the fast dynamical system {x,y) 
which has been studied in Ref. Q. The analysis of the two- variable system 
shows Q that, with T and K fixed, a sub-critical Hopf bifurcation occurs if we 
increase the term H+I (that is, the bias parameter H and the external current 
/). With / = 0, the autonomous behavior is controlled by the bias term H and 
changes dramatically at the bifurcation value = Tatanh(a;^) — {1 — K)xf, 



Close to the bifurcation line, the behavior is that of an excitable cell j5|. Since 
the Hopf bifurcation is sub-critical ("hard excitation"), this model belongs 
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Fig. 1. Different dynamical behaviors of the three-dimensional map for K = 0.6. 
Unless otherwise stated, T = 0.35 and 5 = \ = 0.001: a) cardiac-like spikes (spikes 
with plateau, T = 0.25, xr = —0.5); b) chaotic behavior (T = 0.322, xr = —0.4); 
bursting behavior in c) {xr = —0.45) and d) {xr = —0.6); e) regular spiking 
or "singlet bursting" {5 = X = 0.003, xr = —0.62); f) sub-threshold oscillations 
(T = 0.45, XR = -0.5) and g) fast spiking (T = 0.45, xr = -0.2). For better 
visualization, the curves have been shifted vertically by 2 (dimensionless) units. 
The dotted lines correspond to x = 0. 

to the class of the FitzHugh-Nagumo and Hodgkin-Huxley models (class I 
neurons) [13, 0|, where the stable fixed point is a focus (perturbations lead 
to spiraling behavior in the phase plane). These small oscillations are de- 
scribed in biophysical models, after linearization about the fixed point, as a 
phenomenological inductance behavior (see chapter 10 in [^). 



3 Stability Diagrams for the Three-Dimensional Map 

The three-variable model presents several qualitatively different dynamical 
behaviors which are usual in real biological neurons (Fig. 1). In particular, 
the slow current z{t) is responsible for the bursting dynamics, as discussed by 
Rinzel and Ermentrout . This slow dynamics is controlled by the parameters 
(5, X^xr). In what follows, we focus on the analysis of the stability diagrams 
in terms of these parameters. The lines where the fixed point {xt^^y,,, z^) loses 
stability (for / = 0) have been determined by standard linear stability analysis, 
which yields the equation 

- + [a + (1 - 5)]K^ -a[X + K + l-5]K + Ka{l - 5) = (2) 
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Fig. 2. Left panel: bifurcations for 5 
fixed point is stable below the lines. 



0. Right panel: bifurcations for A = 5. The 



for the eigenvalues A, where a = (1 — x1)/T. 



3.1 The case 6 = Q 



The case 5 = provides a Umiting behavior where the neuron can be "per- 
fectly adaptive". This means that, when the fixed point is stable, adding an 
external current / produces a transient response: firing only occurs while z{t) 
is adapting (phasic firing) 6|. The fixed points can be determined exactly in 
that case: x* = = xr and 2;* = Tatanh(xi?) — (1 — K)xr. This system 
presents monotonic behavior as far as loss of stability is concerned, as shown 
in the left panel of Fig. 2: with all the other parameters frozen, the fixed point 
loses stability by either increasing xr, increasing A or decreasing T. 



3.2 The case \ = 5 



The subspace \ = 5 reflects the simple case where the time scales of the inflow 
(A) and outflow {5) of the z current are the same. As opposed to the 5 = case, 
the bifurcation line in the Xr vs. T plane now presents a minimum, signaling 
that for sufficiently low values of xr, the fixed point can become stable by 
lowering T (Fig. 2d). This non-monotonic behavior in the parameter T is also 
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observed in the planes (T, A) and {X,xr) (Figs. 2e and 2f). 
3.3 Adiabatic approximation 

In the hmit of 5, A ^ 1 one can obtain analytical results by means of an 
adiabatic approximation. By taking 2r to be a quasi-static variable, in the 
fixed point region the stationary value is determined by = {X/6){xr — 
X*). It is thus clear that the parameter xr is essential for controlling the 
equilibrium value of z. Applying the results for the fast sub-system |5|, the 
adiabatic approximation {z^..)^ = He defines a critical value for the rever- 



sion potential: x^(T) = xf + {6/\)Hf = ± [1 - (5/A)(l - K)] ^1 - T/K ± 



(5/A)T atanh 1 — T/fCj . This critical line is plotted in the {xr,T) plane 
for K = 0.6 (Fig. 2d), and is visually indistinguishable from the results for 
X = S = 0.001. 

3.4 The (I,T) plane 

As observed in the previous section, within the adiabatic approximation the 
critical current Ic for the three-variable model can be readily obtained from 
the two-dimensional result, 1^ = — z^ = H^iT) + (A/5)(xc(T) — xr), with 
He = min{H^,H~} and Xc = min{x^,x~}. Remembering that the critical 
x% for / = is Xc + {6/X)Hc, we finally get /^(T) = {X/6){x%{T) - xr). This 
means that the fixed point stability boundary Ic{T) is exactly the same as the 
Xr versus T diagram (Fig. 2d) with a linear shift in the I scale given by the 
above equation. For example, if we set 5 = A ^ 1 and xr = —1, the element 
starts firing at Ic{T) = x^^{T) + 1. 

3.5 The {X,5) plane 

Relaxing the constraint X = 5 and moving away from the adiabatic condition 
A, 5 <^ 1 gives rise to an interesting interplay between 5, A and T. Fig. 3a 
shows the 5 vs A stability boundaries for T = 0.34. While the effect of xr re- 
mains monotonic (i.e. larger values of xr tend to destabilize the fixed points), 
the curves 5c(A) have a maximum. An interesting phenomenon occurs if T is 
lowered (Fig. 3b, with T = 0.3): the shape of the boundaries changes dramati- 
cally, with larger values of A inducing the instability of the fixed point. Fig. 3c 
shows how the curves change as T varies. For larger values of T, the stability 
of the system is weakly dependent on A, depending rather on sufficiently small 
values of 5. The change in the shapes can be intuitively understood on the 
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Fig. 3. Hopf bifurcations in the (A, 6) plane for K = 0.6. The fixed point is stable 
(a) under the lines; (b) to the left of the lines, (c) Evolution of the curves from (a) 
to (b). 

basis of Eqs. 1: large T means small increments in the values of x, rendering 
the 2-equation strongly dependent on 6 (Fig. 3a). For small T, on the other 
hand, S becomes less relevant: since the changes in x are larger, A must be 
small for the fixed point to remain stable (Fig. 3b). 



4 Conclusions 



We studied a three variable nonlinear map which includes a slow dynamics 
over a sub-critical Hopf bifurcation generically found in class I neurons. This 
element is used to represent different behaviors observed in biological cells 
such as excitability, regular spiking, spikes with plateau and bursts [6j. We 
have presented bifurcation lines in different planes in the parameter space. 
This gives us a global understanding of the dynamics of the map and enables 
the proper choice of parameters for locating the excitability region. These maps 
can be coupled by properly modeling chemical or electrical synapses We 
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are currently applying this coupled map framework to large scale biologically 
realistic architectures, results will be published elsewhere. 
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